Hyperbolicity of exact hydrodynamics for three-dimensional linearized Grad's equations 
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We extend a recent proof of hyperbolicity of the exact (to all orders in Knudsen number) linear hydrodynamic 
equations [M. Colangeli et al, Phys. Rev. E (2007)] to the three-dimensional Grad's moment system. A proof 
of an _f/-theorem is also presented. 
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I. INTRODUCTION 

Derivation of hydrodynamics from a microscopic descrip- 
tion is the classical problem of physical kinetics. The 
Chapman-Enskog method [1] derives the solution from the 
Boltzmann equation in a form of a series in powers of Knud- 
sen number e, where e is a ratio between the mean free path of 
a particle and the scale of variations of hydrodynamic fields. 
The Chapman-Enskog solution leads to a formal expansion 
of stress tensor and of heat flux vector in balance equations 
for density, momentum, and energy. Retaining the first or- 
der term (e) in the latter expansions, we come to the Navier- 
Stokes equations, while next-order corrections are known as 
the Burnett (e^) and the super-Burnett (e^) corrections [1]. 

However, as it was first demonstrated by Bobylev for 
Maxwell's molecules [2], even in the simplest case (one- 
dimensional linear deviation from global equilibrium) the 
Burnett and the super-Burnett hydrodynamics violate the ba- 
sic physics behind the Boltzmann equation. Namely, suf- 
ficiently short acoustic waves are increasing with time in- 
stead of decaying. This instability contradicts the i?-theorem, 
since all near-equilibrium perturbations must decay. This cre- 
ates difficulties for an extension of hydrodynamics, as derived 
from a microscopic description, into a highly non-equilibrium 
domain where the Navier-Stokes approximation is inapplica- 
ble. 

Recently, Bobylev suggested a different viewpoint on the 
problem of Burnett's hydrodynamics [3]. Namely, violation 
of hyperbolicity can be seen as a source of instability. We re- 
mind that Boltzmann's and Grad's equations are hyperbolic 
and stable due to corresponding ff-theorems. However, the 
Burnett hydrodynamics is not hyperbolic which leads to no 
_ff -theorem. Bobylev [3] suggested to stipulate hyperboliza- 
tion of Burnett's equations which can also be considered as 
a change of variables. In this way hyperbolically regularized 
Burnett's equations admit the iJ-theorem (in the linear case, 
at least) and stability is restored. 

Inspired by this study, in our recent paper [4] (referred as 



CKK hereafter), we have considered the simplest nontrivial 
example - linearized Grad's moment equation in one spatial 
dimension - and demonstrated that, upon a certain transfor- 
mation, the exact (to all orders in Knudsen number) hydrody- 
namic equations are manifestly hyperbolic and stable. Thus, 
the first complete answer to what is the structure of the ex- 
tended hydrodynamics was obtained. 

In this paper, we extend the CKK result to three- 
dimensional linearized Grad's equations. In addition we prove 
the existence of an i7-function. The paper is organized as 
follows: In Sec. II, through a Dynamic Invariance Principle 
[6-8], we derive equations of linear exact hydrodynamics. In 
Sec. Ill we demonstrate that exact hydrodynamic equations 
are manifestly hyperbolic and dissipative. Then, In Sec. IV 
we stress explicitly how the stability of hydrodynamic equa- 
tions, and therefore the existence of an i/-theorem, arises as 
an interplay between these two basic ingredients of result- 
ing hydrodynamics: dissipativity and hyperbolicity. Finally, 
a conclusion is given in Sec. V. 



II. HYDRODYNAMICS FROM THE LINEARIZED GRAD 
SYSTEM 

A. Linearized Grad's equations in fc-space 

The thirteen moments linear Grad system consists of 13 lin- 
earized PDE's giving the time evolution of the hydrodynamic 
fields (density p, velocity vector field u. Temperature T) and 
of higher order distinguished moments: five components of 
the symmetric traceless stress tensor cr and three components 
of the heat flux q [5]. 

Point of departure is the Fourier transform of the linearized 
three-dimensional Grad's thirteen-moment system: 
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where fe is the wave vector, pk, Uk and Tk are the Fourier com- 
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ponents for density, average velocity and temperature char- 
acterizing deviations from the equilibrium state, respectively, 
and CT/j and q^. are the nonequilibrium traceless symmetric 
stress tensor {'a' = a) and heat flux vector components, re- 
spectively. The overline bar denotes the traceless symmetric 
part of a 2nd rank tensor a, 1i = i (a + a^) — itr(a)I with 
unity matrix I. The system (1) provides the time evolution 
equations for a set of hydrodynamic (locally conserved) fields 
[p, M, T] coupled to the nonhydrodynamic fields cr and q. The 
goal is to reduce the number of equations in (1) and to arrive 
at a closed system for the hydrodynamic fields only. 

To this end, it is common practice to decompose the vec- 
tors and tensors into parallel (longitudinal) and orthogonal 
(lateral) parts with respect to the wave vector, because the 
fields are rotationally symmetric around any chosen direc- 
tion k. We introduce a unit vector in the direction of the 
wave vector, = k/k, k = \k\, and the corresponding 



% Bk + qt, and 
= Q,ek-qi = 0, 



decomposition, Uk 
o-fe = |cr| 'efeefc' + 2crl, where Sk 
and efcgfe : cr^ = 0. 

Upon inserting the above decomposition into (1), and using 
identities, 6^6^ • Cfc = (2/3)efe, Cfeefe : 'ckek = 'ekBk : 
'efcCfe' = 2/3, we obtain the following two closed sets of 
equations for the longitudinal and lateral modes. 
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Equations (2) and (3) are a convenient starting point to de- 
rive closed equations for the hydrodynamic fields. To this end, 
the Chapman-Enskog method amounts to eliminating the time 
derivatives of the stress tensor and of the heat flux in favor of 
spatial derivatives of the hydrodynamic fields of progressively 
higher order. It had already been noted earher [8] that we can 
express the stress tensor and the heat flux vector linearly in 
terms of the locally conserved fields by introducing six, yet 
unknown, scalar functions A{k), . . . , Z{k) for the longitudi- 
nal part: 



ikAul - PBpk - fc^CTfc, 



ikXpk + ikYTk 



k^Zv} 



(4a) 
(4b) 



and, respectively, two functions D{k) and U{k) for the 
transversal component. 



ikD ekwj: 
-k^Uui, 



(5a) 
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where the expressions for the longitudinal components share 
their form with the one-dimensional CKK case. Note that the 
functions introduced should be regarded as exact summation 
of the Chapman-Enskog expansion which amounts to expand- 
ing these functions into powers of fc^ and deriving coefficients 
of this expansions from a recurrent (nonlinear) system, cf. 
CKK and [8]). We do not dwell on this here since we shall 
use a more direct way to evaluate functions A, Z, D,U in 
the sequel. 

Finally, using expressions (4) and (5) in (2), (3) and denot- 
ing as Xk — {pk,ul,Tk,u-j^) the vector of the hydrodynami- 
cal variables, the equations of hydrodynamics can be written 
in a compact form using a block-diagonal matrix Mk, 



dtXk=MkXk, = i Ml ) ^ 



with 

Ml- 

and 



-ik 

-ik{\-k'^B) k'^A 
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Mi=k^Dfi ; 



(7) 
(8) 



where the unit matrix is written in an (arbitrarily) fixed basis 
in the two-dimensional subspace of vectors ti^. As follows 
from an immediate comparison with CKK, and due to the ap- 
parently useful notation, the matrix Ad^ providing the evolu- 
tion of the longitudinal modes, is exactly identical with the 
corresponding matrix (denoted as M in CKK) for the one- 
dimensional case, where lateral modes are absent. The twice 
degenerated transversal (shear) mode is decoupled from the 
longitudinal modes. As a direct consequence, also the in- 
variance equations to be discussed next, which will provide 
us with a set of nonlinear algebraic equations for the un- 
known functions A-Z, divide into two sub-blocks which can 
be solved separately. 



B. Invariance Equations 

In order to evaluate functions A, . . . ,Z,D,U,wc make use 
of the dynamic invariance principle (DIP) [6-8]. Making use 
of DIP in just the same way as for the one-dimensional case 
(CKK) leads to two independent sets of invariance equations 
for the functions A{k)-Z{k). We find that the first set (six 
coupled quadratic equations for A, B, C and X, Y, Z) is iden- 
tical to the one already presented, cf. CKK, Eq. (17). 
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FIG. 1 : Real parts of coefficients Ato Z solving the invariance equa- 
tions, CKK (Eq. 17) supplemented with (10). 



For the transversal modes, the invariance condition reads. 



dot 



{-ikek ■ (ri) = dtqi, 



(9) 



where the time derivative in the left hand side is evaluated by 
chain rule using dtW^. Substituting the functions (5) into (9), 
and requiring that the invariance condition is valid for any wj^, 
we derive two coupled quadratic equations for the functions D 
and U which can be cast into the following form: 



ISfc^D^ + 25PD^ + (10 + 21A;2)£) + 10 = 0, 
3D 



U = - 



2 + 3k^D' 



(10) 



Solution of the cubic equation ( 1 0) with the initial condition 
£'(0) — —1 matches the Navier-Stokes asymptotics and was 
found analytically for all k. This solution is real-valued and 
is in the range D (A:) G [-1.04, 0], whereas C/(fc) £ [0,2.72]. 
The functions corresponding to the longitudinal part of the 
system have been obtained numerically in CKK. Because D 
and U are real- valued, we show in Fig. 1 the real parts for 
all coefficients, while their nonvanishing imaginary parts still 
coincide with those shown in CKK Fig. 4. 

The dispersion relations uj{k) for the five hydrodynamic 
modes are then calculated by inserting these coefficients into 
the roots of characteristic equation det {Mk — uil) = 0, 
where J is a 5 x 5 unit matrix. Analogously, the dispersion 
relations for the remaining non-hydrodynamic modes follow 
from eight (remaining) eigenvalues of (2), (3) with (4), (5). 
All 13 modes are presented in Fig. 2. The resulting hydro- 
dynamic spectrum consist of five modes: the acoustic mode, 
Wac(fc)> represented by two complex-conjugated roots, the 
real-valued thermal (diffusive) mode, (both modes already oc- 
curring in the one-dimensional case) and a twice-degenerated 



FTG. 2: Dispersion relations u){k) for the linearized Grad's sys- 
tem using projected variables, Eqs. (2) and (3). The five hydrody- 
namic modes (diffusive, twice degenerated shear, and two complex- 
conjugated acoustic modes), as well as the eight non-hydrodynamic 
modes are presented as a function of k. While the acoustic mode 
is complex-valued for all k, the remaining modes become complex- 
valued beyond the two visible bifurcation points (at k'^ ~ 0.2175 
and kc ~ 0.303) For k < k'^ the non-hydrodynamic (3D) modes are 
degenerated two and four times, respectively, corresponding to the 
two and four components of and <t^. 



real-valued shear mode (cf. Fig. 2). The occurrence of a real- 
valued shear mode confirms a more general result: in the Un- 
ear regime, the shear mode never undergoes damped oscilla- 
tions. Same as in the one-dimensional case, a critical point 
in the hydrodynamic spectrum occurs at kc ~ 0.303, where 
the thermal mode intersects a non-hydrodynamical branch of 
the original Grad system. Hence, same conclusions hold here: 
for k > kc, the CE method does not recognize any longer 
the resulting diffusive branch as an extension of a hydrody- 
namic branch. Figure 2 further shows the eight (all degener- 
ated) non-hydrodynamic modes, which in opposite to the one- 
dimensional case (offering two non-hydrodynamic modes) 
also exhibit a critical fc at fc^ » 0.2175. 

To summarize, exact hydrodynamics as derived from in- 
variance condition (or, equivalently, by the complete summa- 
tion of the CE expansion as demonstrated in CKK (cf . also 
[8]) extends up to a finite critical value kc, in full agreement 
with the one-dimensional case. No stability violation occurs, 
unlike in the finite-order truncations thereof. Next, we address 
the question about hyperboUcity of exact hydrodynamics in 
the present three-dimensional case. 



III. HYPERBOLICITY OF EXACT HYDRODYNAMICS 

Distinguishing between the real (TZk) and imaginary (Tk) 
parts of matrix Mk (6), we can write the equation of hydro- 
dynamics conveniently as 



dtXk = [R-k - ilk] Xk, 



(11) 
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^ (Re{Ml) \ (iMMl) 0\ 

(12) 

The system (11) is hyperbolic and stable if we can find a 
transformation of the hydrodynamic fields, a;^ =TkXh, where 
Tfc is a real-valued matrix, such that, for the transformed ma- 
trices M'f. = TkMkT^^ it holds 

(i) TV^. = Re(A1'^) and JjJ. = Im(A1'^) are symmetric, and 

(ii) IZ'i^ has non-positive eigenvalues. 

Due to the block-diagonal structure of (6) as well as to the 
fact that CKK has solved the problem of finding a transforma- 
tion with the desired properties for the one-dimensional case, 
the transformation exists also in the three-dimensional case, 
and has the following form: 

where 7^" is explicitly given by CKK Eqs. (25)-(27) in terms 
of k, A-C and X-Z, and 

Thus, the transformation Tfe (13) symmetrizes A^fe and ren- 
ders the exact hydrodynamic equations manifestly hyperbolic. 
Furthermore, the transform 7^ contains only even powers of 
k, because the same is true for the coefficients A-Z. 

The five eigenvalues Ai_5 of 7^^ (or, equally, of TZk), are 

Ai =0, A2 = k^A, A3 = \k'^Y, A4,5 = k'^D. (15) 

From the analysis of the previous section, where we solved 
for coefficients A, D, and Y appearing in (15), cf. Fig. 1, 
it follows that all the eigenvalues Ai_5 are non-positive for 
all k. Note that the matrix TZ'). is diagonal with the diagonal 
elements (15). 



IV. H-THEOREM FOR EXACT HYDRODYNAMICS 

Finally, the hyperbolic structure straightforwardly implies 
an i?-theorem for the exact hydrodynamics (the same holds 
for any lower order approximation, if they are obtained ac- 
cording to the method presented in CKK). Note that, due to 
linearity of the system (1), the choice of a proper iJ-functional 
is not unique. We follow Bobylev [3], and consider an H- 
function - in terms of the transformed hydrodynamic fields - 
defined as: 

H=IJ [P"{r, t) + u'\r, t) + T'\r, t)] d?r. (16) 



Here, hydrodynamic fields x' (r , are defined through inverse 
Fourier transform of the fields x'^. Note that x'{r, t) are real- 
valued because the real-valued transformation 7^ is an even 
function of k,Tk =T-k. Therefore, 

^=IJ [p'kP'-k + < ■ ^'-k + UTLk] d% (17) 
which we abbreviate as -ff = ^ (x'j., x'_^.). Thus, 

dtH = \{{xi,dtx'_^) + {dtxlx'_^)) 

+ li{x'k,-R-'_,x'_,) + (x'_fe,7^;4)). (18) 

Since is an odd function of k,X'_i. = — 1^, terms con- 
taining T' cancel out, and we have, owing to the fact that IZ' 
is even function of k {TZ'_f^ = T^j^), 



dtH = Y^ / \eK^k\'^(fk < 0. (19) 

s=l 

Thus, we have proved the i7-theorem for the exact hydrody- 
namics for k < kc (at k = kc, the eigenvalues A2 and A3 
become complex- valued, as discussed above). 



V. CONCLUSIONS 

In this paper, we have considered derivation of exact hydro- 
dynamics from linearized three-dimensional Grad's system. 
The main finding is that the exact hydrodynamic equations 
(summation of the Chapman-Enskog expansion to all orders) 
are manifestly hyperbolic and stable, thereby extending the 
previous CKK result [4]. To the best of our knowledge, this is 
the first complete answer of the kind. The study supports the 
recent suggestion of Bobylev on the hyperbolic regularization 
of Burnett's approximation. We have also demonstrated, by a 
direct computation, the if-theorem for the quadratic entropy 
function. 
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